Effects due to backscattering and pseudo-gap features in graphene nanoribbons with 

single vacancies 
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We present a systematic study of electron backscattering phenomena during conduction for 
graphene nanoribbons with single- vacancy scatterers and dimensions within the capabilities of mod- 
ern lithographic techniques. Our analysis builds upon an ab initio parameterized semiempirical 
model that breaks electron-hole symmetry and nonequilibrium Green's function methods for the 
calculation of the conductance distribution g. The underlying mechanism is based on wavefunc- 
tion localizations and perturbations that in the case of the first -n — n* plateau can give rise to 
impurity-like pseudogaps with both donor and acceptor characteristics. Confinement and geome- 
try are crucial for the manifestation of such effects. Self-consistent quantum transport calculations 
characterize vacancies as local charging centers that can induce electrostatic inhomogeneities on the 
ribbon topology. 



Graphene is a stable, two-dimensional, sp 2 -bonded car- 
bon allotrope system with exceptional electrical, me- 
chanical and optical properties^. Particularity stems 
from the almost linear dispersion relation in momentum 
space close to the Fermi level, where valence and conduc- 
tion bands meet exactly at the charge neutrality point. 
Provided that a viable method of controllable band-gap 
engineering can be obtained, graphene can also consti- 
tute a valid alternative for post-Si CMOS technology. 
Graphene's peculiar conduction properties can be highly 
compromised by the presence of disorder in terms of lo- 
cal scattering centers (e.g. defects 3 -) or substrate-induced 
interference^. On the other hand, disorder by means 
of chemical functionalization has been proposed to en- 
hance device functionality by controlling the opening of 
exploitable gaps. Dopants can either act directly on the 
electronic structure^, or induce backscattering effects 
within the transport process 8 -. As a consequence, the 
study of conduction in disordered graphene systems ac- 
quires a double significance: on the one hand it becomes 
crucial for performance-related characteristics, while on 
the other it can be used to address band-gap tailoring 
issues. 

Vacancies can make part of crystalline graphene as pro- 
duction faults of the mechanical or the epitaxial growth 
process. Their presence has been experimentally veri- 
fied by transmission electron microscopy^, whereas struc- 
tural, electronic and magnetic properties have been theo- 
retically investigated at a full quantum scaled—. Mod- 
eling within the nearest-neighbor tight-binding formal- 
ism is straightforward: a local point potential U — > oo 
is introduced on the vacancy site that forbids hoppings 
to and from neighboring sitesi -. As a result a semi- 
localized state appears at the Fermi level of the system 
with a C3 point symmetry and a local spi n 10 ' 13 . The po- 
sitioning of this zero-energy mode is highly related to the 
electron-hole symmetrical description of the bandstruc- 
ture under this approach. By considering a next-to-near 



neighbor in the same type of model, valence and conduc- 
tion band mirror-symmetry breaks, and the defect state 
looses its high-symmetry allocation. Accurate ab initio 
calculations show that the exact resonance of this mode 
is located at energies lower than the charge neutrality 
pointii. Pereira et alJ^- have argued that the breaking of 
electron-hole symmetry in the presence of a point defect 
can induce further alterations in the electronic structure 
by means of other semi-localized eigenstates, implying 
that the perturbation induced by the point defect is not 
constrained to a single energy but it expands within the 
energy spectrum. Unlikely, a clear allocation of the role 
of such perturbation with respect to the conduction prop- 
erties of graphene-based systems still lacks. 

Under this perspective, the objective of the present 
work is to extensively investigate the importance of wave- 
function localization induced by defects for the elec- 
trical transport properties of confined graphene. Here 
we present a systematic investigation of conduction and 
charging for single- vacancy damaged armchair and zigzag 
graphene nanoribbons (aGNRs and zGNRs respectively) 
of dimensions within the capabilities of modern litho- 
graphic techniques (width up to 4.6 nm). The basis of 
the transport formalism is a first-principles parameter- 
ized semiempirical Hamiltonian that considers atomic in- 
teractions on the evaluation of distance-dependent over- 
lap integrals, hence, introducing further neighbor inter- 
actions in a natural way. Computational results show 
that there is a clear relationship between the position 
of the vacancy, the resonance of the defect-states and 
the local eigenvector value of the corresponding unper- 
turbated system. Coupling between such features can 
give rise to resonant backscattering phenomena during 
the conduction process that in the case of the first n — tt* 
plateau are associated with the opening of pseudogaps, 
similar to ones obtained for p-type impurities 8 . This 
type of phenomenon emerges for all ribbons within this 
study. Notwithstanding this, a vast complexity charac- 
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terizes the results from case to case, since the eigenvec- 
tor nature of the problem imposes a strict relationship 
with respect to the system geometry in terms of confine- 
ment, position and symmetry. Moreover self-consistent 
Schrodinger/Poisson calculations show that vacancies be- 
come charging centers in a nonequilibrium process and 
induce long-range carrier-density inhomogeneities on the 
ribbon topology. 

We consider the convention of Ref. [HI to classify aG- 
NRs (zGNRs) on the number of the dimcr lines N a 
(zigzag chains N z ) across the ribbon width. Electronic 
Hamiltonians are formed on the basis of the extended 
Hiickel theory^, using a double-^ Slater basis set that 
derives by fitting the bandstructure of bulk graphene ac- 
cording to density functional calculations (see Ref. [H for 
a detailed description of the parameterization). Energy 
eigenstates e Q and their respective eigenfunctions ^ a are 
obtained through a direct diagonalization of the Hamil- 
tonian matrix imposing periodic boundary conditions in 
supercells up to 8.5 nm long. Quantum transport is stud- 
ied within the non-equilibrium Green's function formal- 
ism (NEGF) coupled to the Landauer-Buttiker approach 
for the calculation of the conductance gH. Simulations 
take place within the ballistic regime. Contacts are ideal, 
i.e. of the same width N a (N z ) as the device without the 
presence of vacancies. Edges are passivated with single 
hydrogens. 

We start this study with electronic structure, conduc- 
tance and density of states (DOS) spectra calculated for a 
semimetallic 38-aGNR. Here the vacancy is introduced in 
various positions of a chain transversal to the longitudinal 
axis of the system (see fig. [TJl) . Electronic structure cal- 
culations show that the universal response of the defect is 
to give rise to quasilocalized states^ along nonlocalized 
ones, that due to electron-hole disparity do not preserve 
a mirror symmetry with respect to the charge neutrality 
point. The level of localization for these states can vary 
from strong to weak while their energy resonance highly 
depends on the position of the defect site. When it comes 
to conduction, there are two distinct groups of behaviors 
obtained on the basis of the location of the vacancy sites. 
For sites of the first group, significant conductance dips 
and pseudogap features are present within the valence 
band of the first ir — it* plateau of the system. The 
correspondent DOS interestingly shows a smooth bell- 
like region which expands throughout the conductance 
dip zone. The presence of a non- negligible DOS where 
a conduction gap takes place is strongly correlated with 
resonant backscattering phenomena during the transport 
process. Such feature is strikingly similar to analogous ef- 
fects calculated for p-type (nominally boron) impurities 8 . 
On the other hand, second group sites practically leave 
the pure structure's g — 1 (in units of 2h 2 j e) conductance 
plateau unaltered. Here the first vacancy state below the 
charge neutrality point gives rise to a DOS with a sharp 
peak that decays exponentially and provokes a strongly 
localized perturbation in the electronic structure. A sys- 
tematic analysis for the sum of possible vacant sites shows 
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FIG. 1: (a) Conductance g and DOS as a function of energy 
for a 38-aGNR with single vacancies at different positions N a ■ 
(b) Projection of the wavefunction corresponding to first state 
below the charge neutrality level (upper) and the first state 
above the charge neutrality level (lower) for a non-defected 
38-aGNR. 



that two thirds of the total sites belong to the first group, 
while the remaining one third sites conform with the sec- 
ond group. This categorization strictly stands only for 
the first 7r — 7r* region, since for energies that are more 
distant from the charge neutrality point, further diver- 
gences occur in the conduction properties according to 
purely geometrical criteria. 

The 1/3 — 2/3 relationship presented above can be un- 
derstood with a careful examination of the ribbon's wave- 
functions in the ideal case, i.e. when no alteration of the 
atomic structure is present (see fig. [Tb). The form of 
the eigenvectors within the the entire first conductance 
plateau preserves a 1/3 — 2/3 schema where 2/3 of the 
sites have a finite eigenvector value while the rest 1/3 
have an extremely small value instead. If the vacancy is 
introduced in a position corresponding to the 2/3 group, 
eigenvector symmetry breaks and defect-mode perturba- 
tions spread also to neighboring energies. The sum of 
perturbative behaviors due to the quasilocalized states 
within this plateau gives rise to an increased DOS dis- 
tribution with respect to the ideal case below the charge 
neutrality point. This expanded perturbation is the rea- 
son for electron backscattering during conduction. On 
the other hand, if the vacancy is introduced in the 1/3 
group sites, the perturbation induced remains localized 
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FIG. 2: a) Conductance g as a function of energy for a 
semimetallic 44aGNR with a single vacancy at N a = 22 (left) 
and N a = 21 (right), a) Conductance distribution g as a func- 
tion of energy for a semiconducting 34aGNR with a single 
vacancy at N a = 16 (left) and N a = 17 (right), (c) Projec- 
tion of the wavefunction corresponding to first state below the 
charge neutrality level (upper) and the first state above the 
charge neutrality level (lower) for a non-defected 34-aGNR. 



since the 2/3 symmetry does not break, and no general- 
ized repercussions are inferred in the conductive capacity 
of the system. This picture is not only present in the pre- 
vious aGNR but reflects a general situation for the con- 
duction properties of vacancies in semimetallic nanorib- 
bons with the armchair confinement (e.g. see fig. [5^,). 
A systematic data analysis for aGNRs up to N a — 44 
dimer lines leads to the following empirical rule (which is 
supported by analytical calculations of the first plateau 
wavefunctions for all metallic aGNRsAS) : for N a = 3p + 2 
dimer lines (Vp S N), vacancies at the N a — 3q sites be- 
long to the 1/3 group (Vg G N, < p), while the rest make 
part of the second group. Apart from the semimetal- 
lic ribbons, semiconducting aGNRs behave within a sim- 
ilar framework that infers a less pronounced effect on 
the transport mechanism due to the intrinsic electronic 
bandgap (see fig. [2]d). Unlike though their semimetallic 
counterparts, wavefunctions of semiconducting aGNRs 
do not globally preserve the 1/3 — 2/3 symmetry within 
the first conductance plateau (fig. [2b). Notwithstanding 
this, the mechanism that gives rise to DOS perturbations 
and therefore to conductance dips remains the same, even 
if mainly for vacancies towards the center of the ribbons 
backscattering effects are always present and differences 
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FIG. 3: (a) Conductance g as a function of energy for a 22- 
zGNR with a single vacancy at the N z = 11 (up- left), N z = 10 
(up-right), N z = 7 (down- left) and N z = 6 (down-right) posi- 
tions, (b) Bandstructure of an ideal 22-zGNR with periodicity 
L within the spin-restricted extended Hiickel model. 



are attenuated. 

The case of defected zGNRs is intrinsically more com- 
plicated than that of aGNRs due to the presence of edge 
magnetism (that gives rise to a secondary bandgap 15 ) 
in conjunction with the local spin attributes of va- 
cancies themselves^. On the other hand recent ab 
initio calculations have shown that the most stable 
hydrogen-terminated zig-zag ribbons are not monohydro- 
genated and loose their magnetic ground stated. Here 
we neglect spin interactions and focus on the effect of 
vacancy-induced pseudogaps within the conduction pro- 
cess, which is complementary to but of a higher order 
with respect to the intrinsic bandgap (see Ref. [2(] for a 
discussion on spin-polarized transport). Fig. [3] shows 
representative conductance results for a 22-zGNR with 
single vacancies at four different internal positions of the 
ribbon along a chain transversal to the longitudinal axis 
of the system. The 1/3 — 2/3 symmetry seen in the case 
of aGNRs is not confirmed, reflecting the different form 
of zGNR wavefunctions. On the other hand symmetry 
effects with respect to parity, as well as important dif- 
ferences with respect to the positioning of the various 
defected sites can be obtained. In detail, vacancies at 
the central region of the zGNR show a p-type conduc- 
tion gap that is greatly compromised by the presence of 
the g — 3 plateau area due to band-bending close to 
the Fermi level near the Brillouin zone boundaries (see 
Fig. [HJd). Moving towards the ribbon's edges, even-iV z 
sites become completely metallic whereas odd-A z sites 
give rise to a simultaneous transport gap in both con- 
duction and valence bands 21 . This behavior is similar to 
impurity backscattering with both donor and acceptor 
characteristics, whereas analogous effects have been cal- 
culated for substitutional boron atoms^. The key issue 
for understanding resonant backscattering phenomena in 
zGNRs can be traced back in the high local density of 
states concentration within the n — tt* plateau that is 
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FIG. 4: (a) Schematic representation of the simulated 4- 
terminal system. Front and back gates are isolated by a 
1.9 nm thick SiC>2 layer with a relative dielectric constant 
k — 3.9. Armchair graphene nanoribbons are used as channel 
materials with length I = 8.41nm. (b) Maps of the electro- 
static potential $ (in V) along the ribbon surface for an ideal 
22-aGNR (Vds = 0.4V), a 22-aGNR with a single vacancy at 
N a = 13,1 = 4.38nm (Vds = 0.4V), a 20-aGNR with a single 
vacancy at N a = 4, 1 = 4.28nm (Vds = 0.3V) and a 20-aGNR 
with a single vacancy at 7V a = 6, 1 = 4.28nm (Vds = 0.3V). 



not only limited to edge states, i.e. nonlocalized states 
in an ideal zGNR are present in this energy zone. The 
wavefunctions of these states usually maintain a parity 
symmetry for neighboring sites that can loose balance 
while moving from the center towards the edges of the 
ribbon. According to the exact positioning of the va- 
cancy site the perturbation induced can be stronger or 
weaker while it can affect more than one regions at the 
entire plateau, giving rise to separate resonant backscat- 
tering phenomena during the conduction process. 

For a more realistic description of quantum transport 
and the evaluation of nonequilibrium aspects of con- 
duction, we have performed self-consistent simulations 
within the coupled NEGF-Poisson scheme in order to ev- 
idence the manifestation of charging effects related to 
the presence of vacancies in the atomic lattice. To this 
purpose, we have extended the open-source NanoTCAD 
ViDES cod o 22 ' 23 with extended-Hiickel functionality and 
enhanced matrix operations with optimized numerical 
techniques 24 . We consider the 4-terminal geometry of 
fig. Hi and 20/22aGNRs with length I = 8.41nm as 
channel materials. Single vacancies are introduced at 
various positions of the atomic lattice, whereas the rel- 
atively short length of these aGNRs allows for a short- 
range Coulomb scattering study within the ballistic limit. 
Simulation temperatures are set to T = 300K. Gate 



voltages are fixed to zero (Vds = 0V, Vg 2 s = V) while 
drain-to-source bias is gradually raised from Vds = 0V 
to Vds = 0.5V. Maps of the electrostatic potential $ 
along the ribbon surfaces at the end of the self-consistent 
process can be seen at fig. [4Jd. There are two main charg- 
ing aspects related to the presence of vacancies in these 
samples: on one hand the immediate region around the 
defect acquires a charge concentration that manifests as 
an electrostatic peak on the GNR topology. This fea- 
ture appears for all types of vacancies discussed in the 
previous paragraphs. Directionality of such potential is 
also visible in correspondence to the geometrical orien- 
tation of the vacancy. On the other hand though, an 
important aspect arises in the cases where vacancies are 
associated with impurity-like backscattering issues. Here 
the presence of the defect gives rise to a charge-carrier 
inhomogeneity effect on the GNR topology (e.g. see the 
lower panel of fig. Hfc). Moreover, this process is bias- 
dependent, i.e. by increasing Vds a more intense elec- 
trostatic inhomogeneity pattern can be observed. The 
presence of carrier inhomogeneities on the graphene sur- 
face have been widely attributed to charged impurities, 
and under this perspective important experimental fea- 
tures related to the measured conductivity of graphene 
have been justified 2 —. This study indicates that alterna- 
tively, also vacancies can create similar charging disorder 
whose extend should ideally be calculated for flakes of 
dimensions within the range of the microscalc. 

To conclude, this study has extensively investigated 
conduction and charging properties of armchair and 
zigzag graphene nanoribbons with single vacancies. Fo- 
cus has been put on the underlying backscattering mech- 
anism that proves fundamental for transport-related fea- 
tures like the appearance of pseudogaps within the first 
conductance plateau. A front-end consequence arising 
from this mechanism is that vacancies can behave as p- 
type impurities, while additional donor-like behaviors can 
be observed in the case of zGNRs. It has been argued 
that the positioning of the defect-states within the eigen- 
spectrum in conjunction with geometrical considerations 
that shape the system wavefunctions are the origin of the 
presented phenomena. Moreover self-consistent quan- 
tum transport calculations have evidenced that during 
nonequilibrium, vacancies can induce inhomogeneities in 
the electrostatic topology of the ribbons, in accordance 
with similar effects often attributed to the presence of 
charged impurities. On the basis of such assumptions 
some key points need to be discussed. From an applica- 
tion point of view the association of vacancies with con- 
duction gaps can have a big practical impact on the engi- 
neering of mobility gaps in graphene-based systems, since 
vacancies are either present in the atomic lattice from 
the production stage or can be easily obtained e.g. by 
ion irradiation. Vacancy-concentration has to remain low 
though in order to avoid predominant inelastic electron- 
phonon scattering processes and consequently high re- 
ductions of the electron mobility 3 . A preferential exper- 
imental verification is necessary here. From a method- 
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ological point of view an important issue arises with re- 
spect to the resonance of the defect states. An inaccu- 
rate positioning of these modes in the first conductance 
plateau can lead to a dislocation of the pseudo-gap reso- 
nance, or in some cases, to a complete suppression of such 
effect. In this sense, second or higher neighbor atom- 
istic models, or impurity-like calibrations of the vacancy 
site2£ seem more appropriate for quantum transport cal- 
culations in defected graphene systems. 
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